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I ' Abstract 

■ This paper presents a generalization of our earlier work in [19]. In this paper, the two 

(N 1 concepts, generic regular decomposition (GRD) and regular-decomposition-unstable (RDU) 

variety introduced in [T5] for generic zero-dimensional systems, are extended to the case 
where the parametric systems are not necessarily zero-dimensional. An algorithm is provided 
to compute GRDs and the associated RDU varieties of parametric systems simultaneously 
£ — ■ on the basis of the algorithm for generic zero-dimensional systems proposed in [T^. Then 

the solutions of any parametric system can be represented by the solutions of finitely many 
regular systems and the decomposition is stable at any parameter value in the complement 
of the associated RDU variety of the parameter space. The related definitions and the results 
presented in [19] are also generalized and a further discussion on RDU varieties is given from 
an experimental point of view. The new algorithm has been implemented on the basis of 
CZ5 ' DISCOVERER [28] with Maple 16 and experimented with a number of benchmarks from 

the literature. 

Keywords: parametric polynomial system, regular-decomposition-unstable variety, generic 
regular decomposition 

7-H ■ 

. 

Q ■ 1 Introduction 

^— ^ ' As is well known, solving parametric polynomial system plays a key role in many application 

. fields such as automated geometry theorem deduction, stability analysis of dynamical systems, 

robotics and so on. To solve a parametric system symbolically, a basic idea is to transform the 
system into new systems with special structures or properties so that the solutions of the original 
system can be handled via studying the solutions of the new systems, which is relatively easy. 
. Remarkable examples of such methods are the algorithms for computing comprehensive Grobner 

systems (CGS) and comprehensive Grdbner bases (CGB) [H [H [HI QH QZ1 EB] • The methods 
based on triangular decompositions are another kind of such examples [T] [U [TT] [T3J [501 HU HZ1 

Since Wu's work [27], lots of well-known methods based on triangular decompositions have 
been proposed. An essential concept, "regular chain" (or "normal chain"), and algorithms for 
computing regular chain decomposition have been introduced by Kalkbrener and Yang and 
Zhang |32| independently. For parametric systems, Gao and Chou proposed a method in [7] for 
identifying all parametric values for which a given system has solutions and giving the solutions 
by p— chains without a partition of the parameter space. Wang gave an efficient algorithm for 
computing regular system decomposition |221 1231 125) , which is a generalization of regular chain 
decomposition. The concept of comprehensive triangular decomposition (CTD) introduced by 
Chen et al. in [1] is an analogue of the CGS for solving parametric polynomial systems. 

Two new concepts, generic regular decomposition and regular-decomposition-unstable (RDU) 
variety for generic zero-dimensional systems, are introduced in |19j and an algorithm is proposed 
for computing a generic regular decomposition and the associated RDU variety of a given generic 
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zero-dimensional system simultaneously. The solutions of the given system can be expressed by 
finitely many zero-dimensional regular chains if the parameter value is not on the RDU variety. 
The so-called weakly relatively simplicial decomposition (WRSD) plays a crucial role in the 
algorithm, which is based on the theories of subresultants. 

In this paper, the concepts, generic regular decomposition (GRD) and regular-decomposition- 
unstablc (RDU) variety, introduced in |19j for generic zero-dimensional systems are extended to 
the case where the parametric systems are not necessarily zero-dimensional. An algorithm is pro- 
vided to compute GRDs and the associated RDU varieties of parametric systems simultaneously 
on the base of the algorithm for generic zero-dimensional systems proposed in |19j . Then the 
solutions of any parametric system can be represented by the solutions of finitely many regular 
systems and the decomposition is stable at any parameter value in the complement of the associ- 
ated RDU variety of the parameter space. The new algorithm has been implemented on the base 
of DISCOVERER [55] with Maple 16 and experimented with a number of benchmarks from the 
literature [4j [9l [T2J HHJ [16] . Empirical results are also presented to show the good performance 
of the algorithm. In other words, this paper presents a generalization of our earlier work in |19j . 

First of all, we need to introduce the idea proposed in |19j briefly. For a given generic 
zero-dimensional system P with n variables and d parameters, we considered the parameters as 
"constants" and proposed Algorithm RDUForZD for computing a so-called generic regular decom- 
position T of P in K{U) such that V^jy(P) = UxeTrV ^^ (T\It) , where T is a set of regular 
chains. At the same time, the algorithm would obtain a parametric polynomial such that the 
regular decomposition was stable at any parametric point outside the variety (called RDU vari- 
ety) generated by the parametric polynomial. Roughly speaking, "stable at a parametric point" 
means that the regular decomposition remains after we substitute the point for the parameters in 
P and T (see Definition [5]) . As a result, the original generic zero-dimensional system is "solved" 
except for the case where parameters are on the RDU variety. That is why the decomposition is 
called generic regular decomposition. 

Now we would like to show some new ideas of this paper. If the given system is not generic 
zero-dimensional, to obtain a decomposition with similar properties as in the zero-dimensional 
case, we choose to express the solutions of the system by finitely many regular systems |23j instead 
of regular chains. So we need to generalize the concept "generic regular (chain) decomposition" 
introduced in [TS] into "generic regular (system) decomposition" (sec Definitions [21 El and [5]) . For 
solving a positive dimensional system, a natural idea is to view some variables as parameters and 
call recursively the algorithm for generic zero-dimensional systems proposed in |19j . However, 
to prove the correctness of this procedure, we need to study the properties of characteristic sets 
under specifications carefully (see Lemma [5] and Corollary [2]). Besides, it is worth to notice that 
we have two different interpretations for the results computed by Algorithm ZDToRC and both of 
them play a key role in the proof of the correctness (see Lemma [2]). Finally, we give an algorithm 
which, for any parametric system P, computes a finite set TH of regular systems in K[U] [X] and 
a polynomial B 6 A'[C7], such that 

L V l?(77j( P ) = u [T,H]eTOV W7T (T\ J ff); and 

2. for any a E K d \V u (B), V(P(a)) = U [TjH]eTH V(T(o)\fl'(o)) and [T, if] specializes well at 
a for any [T, H] e TH. 

Please see Algorithm H] in this paper for more details. What's more, for different orderings of 
variables, the efficiency of the algorithm can be different and the RDU varieties can be totally 
different. 

At the end of this section, it is worth to point out that the algorithm provided in this paper 
has a different feature compared to some existing algorithms. The algorithm for computing 
regular system decomposition proposed in |22L 123] uses the so-called variable elimination, which 
computes a main branch at first and then gets the other branches one by one. An incremental 
algorithm, introduced in [H [5J for computing regular chain decomposition, computes a regular 
chain decomposition for some polynomials in the given system at first and then intersects the 
other polynomials with the regular chains one by one. The algorithm proposed in this paper 
makes use of a hierarchical strategy. From an experimental point of view, different strategies are 
suitable for different benchmarks (see Section [4]). 
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The paper is organized as follows. Section [5] provides basic definitions and concepts that 
are needed to understand the main algorithm. Section [3] contains the main algorithm, namely 
Algorithm 01 and some relative subalgorithms. Also we review the description of algorithms 
in our former article. Besides, proofs for these algorithms are presented in this section. Some 
illustrative examples, the empirical data and comparison with previous work along with several 
implementation details are presented in Section^] Scction[S]concludes the paper with a discussion 
on our future work along this direction. 

2 Preliminaries 

The following paragraphs give a brief outline of the vocabulary and tools we will be using 
throughout the paper. All concepts without precise definitions can be found in [5J[57J[2T]. K 
and C stand for the feild of real numbers and the field of complex numbers, respectively. 

Suppose {m, . . . , Ud, Xi, ... , x„} is a set of indetcrminates with a given order u± -< ... -< Ud -> 
x\ < ... ~< x n where {ui, . . . ,Ud} and {xi, . . . , x n } are the sets of parameters and variables, 
respectively. Let U = {tti, . . . , Ud} and X = {xi, . . . , x n }. Suppose K is a field and K denotes 
its algebraic closure. Let K[U] be the ring of polynomials in U with coefficients in K and K(U) 
be the rational function field. A non-empty finite subset P of K[U] [A] is said to be a polynomial 
system or system. If P C A[£/][A]\A[A], it is a parametric polynomial system or parametric 
system. If P C A[A], it is a constant polynomial system or constant system. 

For a non-empty finite subset P GA[[/][A] (A[A]), (P)k[u][x] ((P)k[x)) denotes the ideal 

generated by P in A[£/][A] (A [A]) and yj (P)k[u][x] {yj {^)~g[x]) denotes the radical ideal of 

(P)k[u][x] ((P)k [X ])- For an y F in K \ u ] [ x }\{°} (K[X}\{0}) and for any x G A, if x appears in 
F , F can be regarded as a univariate polynomial in x, namely F = Cqx" 1 + Cix m_1 + . . . + C m 
where Co, C±, .. . , C m are polynomials in A[£/][A\{x}] (A[A\{x}]) and Co ^ 0. Then m is the 
leading degree of F w.r.t. x and is denoted by deg(A, x). Note that if x does not appear in F, 
deg(F, x) = 0. If there exists p (1 < p < n) such that dcg(F, x p ) > and for every i (p < i < n), 
deg(F, x^ = 0, then the class of F is p. If deg(F, x,) = for every i (1 < i < n), then the class of 
F is 0. The class of F in K[U][X}\{0} (A[A]\{0}) is denoted by cls F - If cls F > 0, x clgF is the 
main variable of F and is denoted by mvar(_F). Assume that F = Cgx™ + Cix™ -1 + . . . + C m 
where p = cls F > and Co ^ 0, then Co, denoted by If, is the initial of F and x™, denoted by 
rank(_F), is the rank of F. 

For any P C A'[J7][A], V^(P) denotes the set {(a 1 , . . . , ad+ n ) £ K +n \P(ai, . . . , a d + n ) = 
0,VP £ P}. For any B C K[U], V U (B) denotes the set {(ai, . . . , a d ) G K d \B{a u . . . , a d ) = 
0,VS g B}. And V-^jj(P) denotes the set {(oi,...,a„) G K{Jjf'\P{U, a u . . . , a„) = 0,VP G 
P}. For any F G A[?7][A], the coefficients Bi,...,B t of F in X are polynomials in K[U], 
then Y U (F) denotes V u ({Bi, . . . , B t }). Note that for two finite subsets P and H of K[U][X], 
V 1 ^ !7 y(P\H) denotes the set V^7 !7 r(P)\V^7 !7 r(H). Similarly, we can have V(P\H), V^(P\H) 

and V C/ (P\H). dim(V) denotes the dimension of an affine variety V. 

Triangular Set. A non-empty finite set T = {Ti, T 2 . . . , T r } of polynomials in A[[/][X] 
(A[A]) is a triangular set in A[[/][A] (A[A]) if < cIsti < cls-r 2 < . . . < clsT r - For a triangular 
set T in A[C7][A] (A[A]), I T , mvar(T) and rank(T) denote n TeT I T , {mvar(T)|T G T} and 
{rank(T)|T G T}, respectively. The saturated ideal of a triangular set T in A[[/][A] is defined 
as the set {F G A[[/][A]|It s -F G (T)k\u][X\ f° r some positive integer s} and is denoted by 
sat(T)^j c/ jj X j . Similarly, the saturated ideal of a triangular set T in A[A] is defined as the 
set {F G A"[A]|lT S i ;l G (T)^r x i for some positive integer s} and is denoted by sat(T)^r x ,. 

Suppose F G K[U][X] (A [A]) and T is a triangular set in K[U][X] (A [A]), then F is reduced 
w.r.t. T if deg(A, mvar(Ti)) < deg(Tj, mvar(Ti)) for every i (1 < i < r). A triangular set 
T = {Ti, Ti . . . , T r } in A'[J7][A] (A[A]) is a non- contradictory ascending chain in A[[/][A] 
(A [A]) if Ti is reduced w.r.t. {Ti, . . . ,Tj_i} for every i (2 < i < r). A single-element set 
{F} G K[U] ({F} G A) is a contradictory ascending chain in K[U][X] (A [A]) if F ^ 0. 
Remark that an ascending chain is either a non-contradictory ascending chain or a contradictory 
ascending chain. 
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Successive Pseudo Remainder. For two polynomials F and P in P[P][X] (K[X\) and 
a variable x £ A, the pseudo remainder of P pseudo-divided by P w.r.t. x is denoted by 
prem(P, P, x). Particularly, prem(P, P, mvar(P)) is denoted by prem(P, P). For a polynomial 
F e K[U][X] (K[X}) and a triangular set T = {T u ...,T r } in if[£/][X] (K[X]), the successive 
pseudo remainder |32] of F w.r.t. T is denoted by prem(P, T), namely 

prem(P, T) = prcm(. . . prcm(prem(P, T r ), T r _i), . . . , T±). 

For a finite set P C K[U][X] (if[X]), prem(P,T) denotes the set {prem(P, T) | F e P}. 

Successive Resultant. For two polynomials F and P in P[[/][A] (P[X]) and a variable x G 
X, the resultant [32] of P and P w.r.t x is denoted by res(P, P, x). Particularly, res(P, P, mvar(P)) 
is denoted by res(P, P). For a polynomial F £ K[U][X] (K[X\) and a triangular set T = 
{Ti, ...,T r } in if [f/][X] (^[X]), the successive resultant [32] oi F w.r.t. T is denoted by res(P, T), 
namely 

res(P, T) = res(. . . rcs(res(P, T r ), T r _i), . . . , Ti). 

Regular Chain. A triangular set T = {Ti, . . . ,T r } in P[P][A] (if [X]) is said to be a regular 
chain in K[U)[X] {K[X}), if I Tl _^ and for each i(l<i<r), res(I Tl , {P*-i, . • . ,Ti}) ^0. If T 
is a regular chain in if[C/][_X] (if [X]) and mvar(T) = X, T is a zero- dimensional regular chain. 

Regular SysterrQ. "2 Let T C if [[/][A](if [X]) be a regular chain and P G K[U][X](K[X]). If 
res(P,T) ^ 0, then [T,H] is said to be a regular system in K[U][X](K[X}). 

Proposition 1. S // [T,P] is a regw^ar system in K[U][X], then V-^y(T\P) ^ 0. 

Assignment Homomorphism. For each a — (oi, . . . , ad) € P d , a : if [J7][X] — !> P [X] is a 
homomorphism such that <t> a {F) = F(a,X) for all P £ P[P][A] and we denote 4> a {F) by P(a). 
For a non-empty finite set P C if[Z7][X], P(a) denotes the set {P(a)|P G P} and remark that 
P(a)=0ifP = 0. 

Characteristic Set And Wu's Method. An ascending chain C in if [P][A"] is a characteristic 
set of P in if[P][X] if C C (P)k[u][x] and prem(P, C) = {0}. Theorem Q] below is the so called 
well- ordering principle. 

Theorem 1. E3 y /„ 

exists an algorithm which, for an input non-empty finite subset P C 
if [P][X], outputs either a contradictory ascending chain meaning that ~V K ^ (P) = 9, or a (non- 
contradictory) characteristic set C = {C\, . . . , C{\ such that 

V i?(17)( p ) = V^7y(C\I c ) U uLiV 7?M (P U C U {I C J). 

On the base of Theorem [TJ there exists an algorithm, namely Wu's method, for computing a 
finite sequence of ascending chains Ci, C2, . . . , C m (m > 1) in if [P][X] such that 

(1) Ci, C2, • ■ • , C m is a finite sequence of characteristic sets in if [P][A]; 

(2) If m = 1, V-j^uj(P) = 0. Otherwise, suppose § = {Ci|l < i < m and Cj is a non- 
contradictory ascending chain}, then V-^jy(P) = UcesV ^.^ (C\Ic). 

The set of ascending chains {Ci, C2, . . . , C m } above is said to be a Wu's decomposition or 
characteristic set decomposition of P in if [Ci][X ]. 

For any triangular set T in if [U][X], we denote #(A) - #(T) by d(T, X). 

Definition 1. Let P be a parametric system in K[U] [X] and the set {Ci, C2, . . . , C m } of ascend- 
ing chains be a Wu's decomposition o/P in K[U] [X] . i/d(C,;, X) = for every non- contradictory 
ascending chain Ci, P is said to be a generic zero-dimensional system. Otherwise, P is said to 
be a generic positive-dimensional system. 

Definition 2. Let P be a parametric system in K[U][X] and TH = {[Ti, Pi], . . . , [T s , H s ]} be a 
set of regular systems in K[U][X]. TH is said to be a parametric regular system decomposition 
ofP tn K[U][X], i/V w (P) = Uf =1 V 7?M (TAP I )- 

1 The definition of regular system is the same as that introduced in [4] and different from that proposed in 1221 . 
see more details in [4]. 
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Definition 3. El L e t [T, i2] be a regular system in K[U][X] and a G K . J/T(a) is a regular 
chain in K[X], rank(T(a)) = rank(T) and ves(H(a), T(a)) 7^ 0, i/ien ioe say i/iai £/ie regular 
system [T, U] specializes well at a. 

Definition 4. Let P be a parametric system in K[U][X] and TH = {[Ti, JTi], . . . , [T s ,iT s ]} be 
a parametric regular system decomposition of P in K[U][X]. For any a G R d , i/V(P(a)) = 
U| =1 V(Tj(a)\i?i(a))) and [Tj,.ffj](l < i < s) specializes well at a, then TH is said to 6e stable 
at a. 

Remark that the concept of stable generic regular (chain) decomposition is first introduced 
in Q3]. 

Definition 5. Let TH be a parametric regular system decomposition of a given parametric system 
P in K[U][X]. If there is an affine variety V in K with dim(V) < d such that TH is stable at 
any a G K d \V, then TH is said to be a generic regular system decomposition o/P and V is said 
to be a regular-decomposition-unstable variety (RDU) o/P w.r.t. TH. 

Definition 6. Let T be a zero- dimensional regular chain in K[U] [X] and P G K[U] [X] . Suppose 
H and G are two finite sets of zero- dimensional regular chains in K[U][X]. If 

(1) V^(T U {P}) = U H6 iV^(H), 

(2) V^(T\P) = U GeG V^(G), 

then (H, G) is said to be a weakly relatively simplicial decomposition o/T w.r.t. P in K[U][X]. 

In [19| , we gave an algorithm for computing weakly relatively simplicial decompositions. Here, 
we only present its specification but omit the details. 



Algorithm 1. WRSD 

Input: A zero-dimensional regular chain T = {Ti, . . . , T n } in i4T[L/][X], a polynomial 

P G -ftT[t/][X], variables X = {xi, ...,x n } 
Output: [H, G, B], where 

(1) (H, G) is a weakly relatively simplicial decomposition of T w.r.t. P in X 

(2) B is a polynomial in K[U) such that for any a G K \V U (B), the weakly relatively 
simplicial decomposition (H, G) of T w.r.t. P is stabkj^] at a. 



3 Theory and Algorithm 

We first give some notations. Assume that Alg is a name of an algorithm and pi, .. . ,pt 
is a sequence of inputs of this algorithm. If the output of Alg(pi, . . . ,p t ) is a finite sequence 
Si, ... ,Ss, Si is denoted by Alg(pi, . . . ,Pt)i for any i (1 < i < s) and also said to be the ith 
output of Alg(pi, . . . ,p t ). Given a finite set S = {si, . . . , St} and a map <f> on S, op(5) denotes 
the finite sequence Si, . . . , St and map(s — > <f)(s), S) denotes the set (j)(S) = {(f>(s)\s G S}. 

3.1 Wu's Decomposition Under Specification 

The general idea of Wu's method is presented in Section [3J Some results on Wu's decompo- 
sition under specification are given in this section. 

Definition 7. Let Pi be a parametric system in K[U][X] and § = {Ci,...,C m } be a Wu's 
decomposition 0/P1 in K[U][X]. Suppose C = {C/,i, C;^, . . . , C;^} is a subset of S satisfying 
that 

(1) C;.i is a characteristic set 0/P1. 

(2) Ifk>2, Ci t i (2 < i < k) is a characteristic set of Pi = Pj_i U C/^-i U {Ic;,i_i} where 
Ci.i-i G Ci t i-\. 

2 Please see the definition of stable in 1191 . 
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(3) If k = 1, C; ; i lis a contradictory ascending chain. Otherwise, C; ; j (1 < i < k — 1) is a 
non- contradictory ascending chain and C; ; fe is a contradictory ascending chain. 

Then C is said to be a line o/S and P c = {Pi, . . . ,Pfc} the corresponding systems. 

Lemma 1. Let Pi be a parametric system in K[U][X] and § = {Ci, . . . , C m } be a Wu's decom- 
position of Pi in K[U][X]. Let C = {Ci,i, ■ ■ ■ , Ci.k} be a line of § with corresponding systems 

{Pi, . . . , Pfc}. Then for any a G K \V U (C^k) , there exists a polynomial pi G P; such that 
Pi(a) ^ for any i (1 < i < fc). 

Proof. We prove it by induction on the number k of elements of C. If k=l, it means C contains 
only one element. Since C;^ is a contradictory ascending chain, we can assume that C; : i = {C\i} 
where C/^i G K[U] and we know that C/^i G (Pi)- Suppose Pi = {/i,i, . . . , fi,tt}- Then C/,i 
can be written as C;.i = ^2,j—\hjf\ t j where hj G A[[/][X] for any j (1 < j < ti). Thus 

C;,i( a ) = Y^j=i ^i( a )/ij( a ) holds. Since Ci.\(a) ^ 0, there must exist some fi <ei G Pi such that 
/i !Gl (a) ^ 0. Let pi = /i jEl and we are done. 

Now we assume that the conclusion holds when k < N (N > 1). Suppose k = N. 
£2 = {Cj 21 • • • , Cjjfe} C C is a line of Wu's decomposition of P2. According to the induc- 
tion hypothesis, for any a G K d \V u (Q^), there exists a polynomial G Pi such that Pi(a) ^ 
for any i (2 < i < k). As is known to us, P 2 = Pi U U {lc, j}, if Pi G Pi, let pi = j»2 and 
the conclusion holds obviously. Otherwise, P2 G {Ic, 1} U C;^. If P2 G C^i C (Pi), supposing 
Pi = {fi.i, . . . , /i,ti}j then P2 can be written as P2 = S*=i hjfi.j where hj G A"[[/][A] for any 
j (1 < j < ii)- Thus ^2(0) = S*=i ^j( a )/i,j( a ) holds. Since ^2(0) ^ 0, there must exist some 
/i, ei G Pi such that fi iei (a) ^ 0. Let p\ = f\_ &1 and we are done. If P2 = {Ici 1}, it implies 
Ci t i(a) ^ 0. Since C;,i can be written as C;.i = J2j=i9jfi,j where gj G A'[J7][A] for any j 
(1 < j < ti). Thus there must exist some f\ t e t € Pi such that /i. ei (a) ^ 0. Let pi = fi,e t and 
the conclusion holds. 

□ 

With Lemma [I] Corollary [T] holds obviously. 

Corollary 1. Let P be a parametric system in K[U][X] and 8 = {Ci,...,C m } be a Wu's 
decomposition of P in K[U][X\. Let C = {C;^, C;^, • • • , C;^} be a line of S. If P = {p}, 

containing only one polynomial in K[U][X], p{a) ^ for any a G K \V u (Ci t k). 

Lemma 2. Suppose C = {C\, . . . , Ct} is a non- contradictory ascending chain and a character- 
istic set of parametric system P in K[U\[X\. For any a G K , V(P(a)) = V(C(a)\Ic(a)) U 
Uj =1 V((P(o)UC(o)U{I Cj (a)})). 

Proof. For any a G K , if Ic(a) = 0, the conclusion holds obviously since Ic(a) = Ili=i ^Ci( a )- 
Now we prove the conclusion when Ic(a) ^ 0. According to the definition of characteristic 
set, wc know that C C (P) and for any p G P, Ici fcl ■ • Tc t fc 'p = <?iCi + • • • + q t C t where 
q l G K[U][X] for any i (1 < i < t). Then C(a) C (P(a)) and l Cl (a) kl ■ ■ ■ l Ct {a) kt p(a) = 
qi (a)Ci(a) + ■■■ + q t (a)C t (a). Therefore V(C(o)\I c (a)) G V(P(a)) G V(C(o)) and V(P(o)) = 
V(C(o)\I c (o))UV(P(o)U{I c (o)}). Since C(a) C (P(a)), V(P(a)U{I c (a)}) = V(P(a)UC(a)U 
I c (a)) = U* =1 V((P(a) U C(a) U {I Ci (a)})). We arc done. □ 

According to Theorem [T] and Lemma we can get Corollary [3] easily. 

Corollary 2. Let P be a parametric system in K[U][X] and {Ci, . . . , C m } be a Wu's decom- 
position ofP in K[U][X]. Suppose S = {C;|l < i <m and C,; is a non- contradictory ascending 
chain} and C§ = {Cj|l < i < m and Ci is a contradictory ascending chain}. Then for any 

a G K d \(U CSeCS V u (CS)), V(P(o)) = U Ce sV(C(a)\Io(a)). 

3.2 Converting To Regular Systems 

Let T be a triangular set in A[[/][A]. Wc can compute a set of finite sequence of regular 
systems TH = {[Ti, H{\, . . . , [T s , H s }} in K[U][X] and a polynomial B G K[U] on the basis 
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of Algorithm [U such that V-^ 7 j(T\I T ) = U| =1 V^7^-(T<\i^) and for any a G K d \X u {B), TH 
specializes well at a and V(T(a)\Ix(a)) = Uf =1 V(Tj(a)\i?i(a)). The algorithm is presented as 
Algorithm EH which plays a key role in Algorithm 2] proposed in the next section. 

Algorithm [5] below was proposed in [T!5] for zero-dimensional case. We just give its specifi- 
cation here. Our focus in this paper is how to deal with the case where the triangular set is 
positive-dimensional. So, we need to convert a triangular set with mvar(T) C X to a set of 
regular systems. 



Algorithm 2. ZDToRC 
Input: A triangular set T in A[[/][A"] with mvar(T) = X, variables X = {x±, . . . ,x n }. 
Output: [G, B], where 

(1) G is a finite set of zero-dimensional regular chains in A[[/][A] such that 
V W (T\I T ) = U GeG V w (G); 

(2) B is a polynomial in K[U) such that for any a G K \Y U (B), 
V(T(a)\I>r(a)) = U GeG V(G(a)) and G specializes well at a for any G € G. 



According to Algorithm^ the following Proposition [5] is clear. 

Proposition 2. Let T be a triangular set in K[U][X] and ZDToRC(T, mvar(T)) = [G,B]. If 
G ^ 0, I T (o) + for a <£ X U (B). 

Now suppose T = {Ti,...,T/} is a triangular set in A"[[/][A] and mvar(T) C X where 
variables X = {x\, . . . ,x n } and parameters U = {iti, . . . , ltd}. It is interesting to show that 
we have two versions to interpret the relationship between T and the results computed by 
ZDToRC(T, mvar(T)). Assume that ZDToRC(T, mvar(T)) = [G, B]. Let £>(T) = mvar(T) and 
J"(T) = A\mvar(T). On one hand, T can be regarded as a triangular set in K[U, J"(T)][S(T)]. 
At this point, according to Algorithm [2J we know that 

( 1 ) V X((7,B(T))( T \ I t) = U G GGV x(;7iJr(T)) (G), 

(2) for any a e K^^XV^^iB), V(T(o)\I T (o)) = U GeG V(G(a)) and G specializes well at 
a for any G G G if G ^ 0. 

On the other hand, T can also be regarded as a triangular set in A[J7][J r (T)][B(T)]. Let 
JC = K{U), U = J"(T) and X = B(T). Then according to Algorithm^ 

(3) V W (T\I T ) = U GeG V^(G), 

(4) for any a G lC n ^ l \V u (B), V^(T(a)\I T (a)) = U GeG V^(G(a)) and G specializes well at a for 
any G G G if G ^ 0. 

Remark that the above statements (1) and (3) are exactly the same since K(U,U) = K{U). As 
discussed above, we have the following lemma by statements (2) and (4) 

Lemma 3. Suppose T is a triangular set in K[U][X] and ZDToRC(T, mvar(T)) = [G,-B]. Then 
V^(T\I T • B) = U GeG V 7?M (G\B) and V(T(o)\I T (o) ■ B(o)) = U GeG V(G(a)\B(a)) for any 

aeK d \Y u (B). 

Theorem 2. Algorithm^ terminates correctly. 

Proof. For a given triangular set T = {T 1 , . . . , T t } in K[U] [X]. Suppose ZDToRC(T, mvar(T)) = 
[G , Bo]- 

Firstly, we prove Algorithm^ terminates. If Bq G A^[C7], it terminates obviously. Otherwise, 
B e K[U][jF{T)] c K[U][X}. Assume that {Ci, C 2 , C m } is the Wu's decomposition of {B } 
in A[J7][X] computed by Wu's method and § = {G|l < i < m and Cj is a non-contradictory 
ascending chain}. Obviously m > 1 since Bq g K[U][J 7 (T)]\K[U}. For any Ci G §, let T^ = 
TUC,. Then we know that mvar(T) C mvar(Tj), which means d(T i: X) < d(T, X). Therefore, 
it is clearly that Algorithm [3] terminates within finite steps of recursion. 

Now we prove the correctness by induction on the recursive depth h. If h = 1, according 
to Algorithm [3l Bq G K[U] and the conclusion follows from Lemma [3l Assume that the con- 
clusion holds for h < N (N > 1). When h = N, Bq G K\U\[F{T)]\K\U). Wc can assume 
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Algorithm 3. TSToRS 



Input: A triangular set T = {T\, . . . ,7]} c if [?7][A], variables X = {xi, . . . ,x n }. 
Output: [G,£>], where 

(1) G is a finite set of regular systems in iT[£/][X] such that 
V^(T\I T ) = U [Gifl]gG V^(G\ff); 

(2) B is a polynomial in K[U], for any a G K \V U (B), 

V(T(a)\Ix(a)) = U[ G £f] 6G V(G(a)\ii(a)) and [G,H] specializes well at a for any 
[G, H] e G if G ^ 0. ' 

1 VF:=ZDToRC(T,mvar(T)) 

2 G:=map(i -> [t, W 2 ], W x ) 

3 if W 2 G if [?7] then 

4 |_ return [G, W 2 ] 

5 Compute a Wu's decomposition {Ci, C2, . . . , C m } of in if[[7][X] by Wu's method 

6 B:=l 

7 for i = 1 —} m do 

8 if Ci is a contradictory ascending chain then 

9 B:=B ■ op(Ci) 

10 else 

11 Tj :=TUC, 

12 G := G U TSToRS(T 4 , X)\, B := B ■ TSToRS(T 4 , X) 2 

13 return [G, B] 



that {Ci, C2, • ■ • , C,„} is the Wu's decomposition of {£>o} in if \U][X] computed by Line 5 in 
Algorithm [3] Let M = {i\l < i < m and Ci is a non-contradictory ascending chain} and 
CM = {1,2,..., m}\M. For any i G M, let T< = T U Q and then T 4 is a triangular set in 
K[U}[X}. Suppose Cj = {C-j} (j G CM), TSToRS(Ti,mvax(Ti)) = [G.,,5,] for any i (i G A/"). 
Then according to Algorithm^ G = U GeGo {[G, B ]} U U ieA ^G, and 5 = H cA -<", ■ ILejV^- 
By Lemma [21 Wu's method and the induction hypothesis, we get 

V W (T\I T ) 

= V W (T\I T • So) U V^y(T U {B }\It) 

- u GeGo v^ I7T (G\i?o) u (v to (t\i t ) n v 7?M (i?o)) 

= u G6 GoV^y(G\Bo) U (V OT (T\I T ) n UeArV^yCCAlcJ) 

= U GeGo V 7?I!7 y(G\Bo) U (U ie ^V w (TAl Ti )) 

= U GeGo V^(G\.Bo) U (U ieA f U [G)H]eGi V^y(G\F))). 

Therefore, the statement (1) in the specification of Algorithm [3] holds. 

For any a e K d \V u (B), Cj(a) ^ and B,(a) ^ for any i £ M and j G CA/". Thus by 
Corollary [21 V(I?o(a)) = Uie/v'V(C,(a)\Ic i (a))- By Lemma [3] and the induction hypothesis, we 
get 

V(T(o)\I T (a)) 

= V(T(o)\I T (o) • B (a)) U V(T(a) U {B (a)}\I T (a)) 

= (U GeGo V(G(a)VB (a))) U (V(T(o)\I T (o)) n V(B (o))) 

= (U GeGo V(G(a)\B (a))) U (V(T(a)\I T (a)) D U l& vV(G( a )\I Cl (a))) 

= (U GeGo V(G(a)\Bo(a))) U (U; eA rV(T;(a)\I Ti (a))) 

= (U G6Go V(G(a)\Bo(a))) U (U l&v U [Gi?]eGi V(G(a)\fT(a))). 

In addition, by Lemma[T]and Corollary[I] Bo(a) 7^ 0. Thus [G,i?o] specializes well at a for every 
G 6 Go according to Algorithm^ By the induction hypothesis, we know that [G,H] specializes 
well at a for every [G, H] G Gi for any i G M. Therefore, the statement (2) in the specification 
of Algorithm [3] holds. □ 
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Remark 1. By Algorithm^ for any regular system [Ti,Hi] in the first output of Algorithm^ 
res(Lrj,Ti) is a factor of Hi. If B 6 K[U] is the second output of Algorithm^ by Corollary^ 
we know that V 1 ^ U j{Tes(l Tll r T l )) C Yj^jjiH,) and V^res^, T 4 )) C V u (Hi) C V U (B). 

3.3 Computing RDU 

We present the main result in this section. Algorithm 0] shows how to compute a generic 
regular system decomposition and the associated RDU of a given system simultaneously. 



Algorithm 4. RDU 



Input: A parametric system P in A[[/][A], variables X = {x\, . . . ,x n }- 
Output: [TH, B], where 

(1) TH is a finite set of regular systems in A [I/] [A] such that 
V i?(!7)( P ) = U [Tiir]GTH V w (T\if), 

(2) B a polynomial in K[U] for any a e K d \V u (B), V(P(a)) = U [T ^ ]eTH V(T(a)\ff(a)) 
and [T, H] specializes well at a for any [T, H] G TH. 

1 Compute a Wu's decomposition {Ci, . . . , C m } of P in A[J7][A] by Wu's method 

2 B:=l 

3 TH:=0 

4 for i = 1 — > m do 
if C; is a contradictory ascending chain then 

| B:=B-op(Ci) 

7 else 

8 W:=TSToRS(Ci, A) 

9 TH:=TH U Wx, B:=B ■ W 2 



10 return [TH, B] 



Theorem 3. Algorithm^ terminates correctly. 

Proof. The termination follows from the termination of Algorithm [3] We only need to show the 
correctness. In fact, the statement (1) in the specification of Algorithm [5] follows from Wu's 
method and Algorithm [3] and the statement (2) follows from Corollary [5] and Algorithm [3J □ 

Corollary 3. Let P be a parametric polynomial system, RDU(P,A) = [TH, B\. Then TH is 
stable at any a <G K \V U (B). 

The proof of Corollary |3j is similar to that of generic zero-dimensional case. For more details, 
please see [15] . 



4 Examples and Implementation 

In this section, we show by examples how our algorithms work. In addition, we run some 
benchmarks and report comparison to other tools with similar function. 

Example [T] is designed to illustrate how Algorithm 0] works. 
Example 1. Consider the system 




(ux + l)z 3 + (vy + l)z 2 + wxz + 1 
ux + I 



where x, y and z are variables (x -< y -< z) and u, v and w are parameters. 

Step 1: According to the first step of AlgorithmHJ we get Wu's decomposition S = {Ci, C2, C 3 } 
of P in IR[u, v, w][x, y, z] where Ci = {ux + 1, u + uvyz 2 + uz 2 — wz}, C2 = {ux+ 1, vy + 1, u — wz} 
and C3 = {-v 2 u 5 w}. 
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Step 2: Let TH = and B = 1. 

Step 3: Because Ci and C 2 are both non-contradictory ascending chains and C3 is a con- 
tradictory ascending chain, we need to execute TSToRS(Ci, [x, y, z]}) and TSToRS(C2, [x, y, z}). 

Step 3.1: According to TSToRS, we execute ZDToRC(Ci, mvar(Ci)) where mvar(Ci) = 
{x, z}. It returns W = [{{ux + 1, u + uvyz 2 + uz 2 u — wz}}, u(vy + 1)]. Since W2 ^ K[u, v, w], 
we execute wasolve(u(vy + 1)) and it returns Cn = {vy + 1}, Ci 2 = {u}. 

Step 3.2: Let Ti = TU {yv + 1}. Now we need to execute TSToRS(Ti). When 
we execute ZDToRC(Ti, mvar(Ti)), it returns [0,-yu]. Since vu 6 K[u,v,w], the output of 
TSToRS(Ci, {x, y, z}}) is [{[{mx + 1, u + uvyz 2 + z 2 ■ u — wz}, u(yy + 1)]}, uvw]. 

Step 3.3: Let TH := THU{[{ua;+l, u + uvyz 2 + z 2 -u-wz}, u(vy + l)]} and B = B-uvw. 

Step 3.4: According to TSToRS, we execute ZDToRC(C 2 , mvar(C 2 )) where mvar(C 2 ) = 
{x, y, z}. It returns W = [{ux + l,vy + l,u — wz}, -tiw]. Since W 2 G K[u, v, w], the output of 
TSToRS(C2, [x, y, z]}) is [{[{wx + 1, vy + 1, u — wz}, —uvw]}, uvw] 

Step 3.5: Let TH := TH U {[{ux + l,vy + l,u- wz}, -uvw]} and B = B ■ uvw. 

Step 3.6: Since C3 is a contradictory ascending chain, we execute B := B ■ —v 2 u 5 w. 

Step 4: Finally, we get B = u 2 v 2 w 2 and T = {[{w.t + l,u + uvyz 2 + uz 2 — wz},u{vy + 
1)], [{ux + 1, vy + 1, u — wz}, —uvw]}. 

In the example, we factor the polynomials and let the polynomial be squarefree in some steps. 

With the result above, we get a regular system decomposition of P in R[it, v, w][x, y, z] and 
a RDU variety V = {(u,v,w) G C 3 \uvw — 0}. According to the specification of Algorithm^ 
for any a g C 3 \V, V(P(a)) = U[ T , ff ] eTH V(T(a)\i?(a)) where [T(a),H(a)] is a regular system. 
With the definition of RDU variety, we get that V is the RDU variety of P w.r.t TH. 

Example [5] is provided by Changbo Chen, which is a good example to show how the ordcrings 
of variables affect the results and the efficiency of Algorithm |U 

Example 2. Consider the parametric system 

d 4 d 3 r + rf - d 4 d 3 r 2 + d\d\ - d 4 d\ - d|rf 3 + d 4 d 3 + Z - r)t 4 + (-2r 2 d 4 r + 2r 2 d\+ 
2r 2 d 4 d\ - 4r 2 d 3 dj + 2r^d 4 + 2r 2 d 4 )t 3 - (2r| - 2r + kd\r\ + Ad 2 4 + 2Z - 2d 2 4 dl)t 2 + 
(-2r 2 d 4 r + 2r 2 d 4 d\ + 2r 2 d 4 + 2r 2 d\ + 4r 2 d 3 dj + 2r\d 4 )t + r\ + d\d 3 - d 4 d 3 r+ 
d 4 d 3 r 2 + Z — r — d 4 d 3 + d\d\ + d 4 d^ 

where r, Z , t are variables and r 2 , d 3 , d 4 are parameters. 

The 3 variables can be ordered in 6 different ways. We tried all these orders when calling 
Algorithm 0] for this example. We are interested in the RDU varieties output by the algorithm, 
so we only report the second output of the algorithm. 

By calling RDU(P, [r, t, Z\) and RDU(P, [t, r, Z]) , we get 

RDU(P, [r,t,Z]) 2 = r 2 d 4 
RDU(P, [t,r, Z]) 2 = r 2 d 4 

By calling RDU(P, [Z, r,t]), we get 

RDU(P, [Z,r,t]) 2 = (-d 3 + d 4 rl + d 4 )d 4 r 2 (d 4 d 3 - 1) 
By calling RDU(P, [r, Z,t\), we get 

RDU(P, [r, Z, i]) 2 = r 2 d 4 (d 3 - d 4 r\ - d 4 ) 

By calling RDU(P, [t, Z,r\), we get 

RDU(P, [t, Z, r]) 2 - r 2 d 4 {d\ - 2>r 2 2 ){d 4 d 3 - 1) 

RDU(P, [Z, t, r]) 2 is too huge to be listed here. It has 11 factors and contains 10838 terms 
(after expanding). 
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Table 1. Timings (in second) of Example^ under different orders of variables 





[r,t,Z] 


[t,r,Z] 


[Z,r,t] 


[r,Z,t] 


[t,Z,r] 


[Z,t,r] 


time 


0.016 


0.015 


0.016 


0.046 


0.016 


3.931 



The timings for the above computation are shown in Table [T] How to choose a suitable order 
in advance is an interesting topic for our future work. 



Table 2. Comparing RDU and Triangularize 



7 

number 


system 


T J 

u 


X 


wusolve 


TSToRS 


total 


Tr i angul ar i z e 


1. 


Hcrcman-2 


1 


7 


0.063 


0. 


0. 063 


0.452 


2. 


Hcrcman-8-8 


3 


5 


0.281 


0.015 


0.296 


0.188 


3. 


Maclane 


3 


7 


0.093 


0.078 


0.171 


0.156 


4- 


MontcsS7 


1 


3 


0.967 


0.015 


0.982 


0.063 


5. 


MontcsSll 


3 


3 


0.016 


0. 


0.016 


0. 


6. 


MontesS12 


2 


6 


0.109 


0.078 


0.187 


0.093 


7. 


MontesS13 


3 


2 


0.016 





0.016 


0.031 


8. 


MontesS14 


1 


4 


0.031 


0.016 


0.047 


0.141 


9. 


MontesS15 


4 


8 


0. 


0. 


0. 


0.016 


10. 


MontesS16 


3 


12 


0.016 


0. 


0.016 


0.156 


11. 


MontesS18 


2 


3 


0.296 


0.031 


0.327 


0.811 


12. 


AlkashiSinus 


3 


6 


0. 


0. 


0. 


0.047 


13. 


Bronstcin 


2 


2 


0.016 


0. 


0.016 


0.062 


14. 


Cheaters- homotopy-easy 


4 


3 


0.358 


10.749 


11.107 


0.047 


15. 


Cheaters-homotopy-hard 


5 


2 


0. 


32.867 


32.867 


0.062 


16. 


Gcrdt 


3 


4 


0.015 


0. 


0. 015 


0.016 


17. 


Lanconelli 


7 


4 


0.047 


0. 


0.047 


0. 


18 . 


Lazard-ascm2001 


3 


4 


0.811 


0.047 


0.858 


0.39 


19. 


Ley kin- 1 


4 


4 


0.312 


0. 


0.312 


0.531 


20. 


Neural 


1 


3 


0.047 


0.016 


0.063 


0.093 


21. 


Pavelle 


4 


4 


0.109 


0.062 


0.171 


0.203 


22. 


SY14 


2 


2 


0.016 


0. 


0.016 


0. 


23. 


Wang93 


2 


3 


0.031 


0. 


0.031 


0.046 


24- 


zhou3 


6 


11 


0.062 


0.047 


0.109 


0.281 


25. 


zhou4 


4 


7 


0.016 


0.016 


0.032 


0.093 


26. 


KdV 


15 


11 


0.406 


0. 


0.40 6 


0.031 


27. 


P3P 


5 


2 





0.031 


0.031 


0.047 


28. 


SBCD23 


1 


3 


0.031 


0. 


0.031 


0.047 


29. 


SBCD24 


1 


4 


0.655 


0.016 


0.671 


0.312 



We have implemented our algorithms with Maple 16. More specifically, Wu's method for 
computing parametric triangular decompositions introduced in Section [5] is implemented as a 
function wusolve and Algorithm [T] is implemented as a function WRSD. We ran many examples 
collected from other papers |15[ 1121 0]. At the same time, we compare the running time with 
Triangulariz^l in RegularChains which can compute regular decompositions of given polyno- 
mial systems with parameters. Throughout this section, all the results are obtained in Maple 
16 using an Intel(R) Core(TM) i5 processor (3.20GHz CPU and 2.5 GB total memory) and 
Windows 7 (32 bit). The empirical data about timings is presented in Table [5J 

In Table [21 the column marked U and X mean the number of parameters and variables, 
respectively. The column marked wusolve means the time used by wusolve at the first step 
in RDU (see Algorithm [4]) . The column marked TSToRS means the time used by TSToRS. The 
column marked Triangularize means the time used by Triangularize. Some data which shows 
0. means that the data is less than 0.001 and ignored by system. 

From Table [5J wc find that the cost of wusolve takes up a majority of the total time in most 
examples when solving practical problems as shown. Comparing with Triangularize, for some 
examples, our algorithm is faster than Triangularize. As for the example Cheaters-homotopy- 
easy and Cheaters-homotopy-hard, the time used by WRSD seems too much. This is not reasonable 
and optimization should be done in the future. 

3 Plcase find more details on Triangularize from the help document of Maple 16. 
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5 Conclusions 



The focus of this paper is how to decompose a parametric system into regular systems at 
"parameter level" and the property of this decomposition under specification. We provide an 
algorithm for computing GRDs of parametric systems and the related RDU varieties simulta- 
neously no matter the systems are generic zero-dimensional or positive-dimensional, which is a 
generalization of our earlier work in |19j for generic zero-dimensional case. Then any parametric 
system in A'[£/][A] can be decomposed into finitely many regular systems and the decomposi- 
tion is stable at any parameter value in the complement of the associated RDU variety of the 
parameter space. That is to say, once we obtain such decomposition, all the solutions of the 
original system are expressed by some regular systems, except for some possible solutions over 
the parameter values on the associate RDU variety. 

Note that, using Algorithm U of generic regular decomposition, we may get a complete de- 
composition of a given parametric system P step by step. A rough procedure is as follows. 
Firstly, for certain variables X and parameters U, we call RDU(P,A) and get a set of regular 
systems and a parametric polynomial B G K[U]. Then, let U\ = U\{u} and X\ = X U {u} 
where u G U and we call RDU(P U {B}, X\) to get a set of regular systems and a parametric 
polynomial B\ G -fsT[C/i]. Continuing this iteration, until \Ji = % and Bi G K for some i, we can 
get a complete decomposition finally. This is what we called hierarchical strategy in Section [1] 
What can we benefit from this strategy? According to this method, we can stop at any step 
of the iteration especially when the computation is hard to be finished. Then we get a generic 
regular decomposition and a set of parametric polynomials which can determinate a RDU variety 
(low dimensional variety). For some huge problems, limited by the computational capacity of 
micro computer, we can get a partial solution, which is useful if one cannot get any information 
by other complete methods. Of course, the procedure should be described clearly and proved to 
be correct. That is one of our future work. 
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